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Galaxy clusters are the most massive systems in the known universe. They host relativistic cosmic ray populations and are thought 
to be gravitationally bound by large amounts of Dark Matter, which under the right conditions could yield a detectable 7-ray flux. 
Prior to the launch of the Fermi satellite, predictions were optimistic that Galaxy clusters would be established as 7-ray-bright 
objects by observations through its prime instrument, the Large Area Telescope (LAT). Yet, despite numerous efforts, even a 
single cluster detection is still pending. 


1. Introduction 


implications for future searches with the LAT. 


Galaxy clusters (GC) represent the largest virialized ob¬ 
jects that are believed to have formed through a hierarchi¬ 
cal build up of structures over the evolution of the uni¬ 
verse. In this picture, baryonic matter accretes towards the 
gravitational well caused by large amounts of Dark Mat¬ 
ter (DM) which make up 26% of the energy density of the 
Universe |Ade et fl/.|2014[ . Through gravitational interac¬ 
tion, smaller structures merge with one another, forming 
groups of galaxies and eventually clusters. Determining 
the exact nature of DM is one of the greatest challenges of 
modern physics and weakly interactive massive particles 
(WIMPs) prove to be strong candidates fulfilling the role 
as DM particle |Bergstrom||1999| [2Q09| |. The neutralino, 
which in several extensions of the standard model of parti¬ 
cle physics is predicted to be the lightest stable supersym¬ 
metric particle, provides a natural WIMP candidate. In 
many of these models the neutralino can self-decay or an¬ 
nihilate into lighter standard model particles, among others 
high energy 7 rays which, if observed, can be used to trace 
back the origin of the interaction and indirectly detect DM 
|Baltz et al. 12008) . 

While clusters are promising targets due to their large 
DM content, predicted 7 -ray emission on top of that of in¬ 
dividual cluster member galaxies constitutes an irreducible 
foreground. This foreground emission arises from cos¬ 
mic ray (CR) interactions with the intra-cluster medium 
(ICM) and is motivated by conventional astrophysics [see, 
e.g. Petr osian e t al. 2008 ( for a review], while the observa¬ 
tion of DM-induced 7 rays may be regarded as a somewhat 
more exotic signature [ Pinzke et <3/. [20111 . 

In this contribution I will review the most recent stud¬ 
ies of GCs with the Fermi Large Area Telescope (LAT) 
undertaken by the instrument team. I will start by summa¬ 
rizing recent efforts aiming at the astrophysical emission 
scenario of CR interactions resulting in a detectable 7 -ray 
flux (Section [2]) and briefly report on work in progress in 
regards to DM constraints that can be obtained from GCs 
(Section [3]). Finally, I will discuss one of the challenges 
involved when searching for large extended sources such 
as GCs (Section [4]) and conclude by commenting on the 


2. Cosmic Ray Induced 7 rays 


The majority of the baryonic mass in GCs is present in 
the form of hot ionized gas, the ICM, which has been de¬ 
tected via thermal X-ray emission observed by contempo¬ 
rary space telescopes such as ROSAT or XMM-Newton 
[see, e.g. Kaas tra et al. ||2008] for a review]. In addition, 
large scale radio synchrotron emission has been detected in 
a number of the most nearby clusters which can be classi¬ 
fied into halos and relics (Ferrari et al. 2008| . The latter ap¬ 
pear polarized, while the former are not, suggesting a dif¬ 
ferent emission mechanism to be at play. The observation 
of radio-synchrotron emission indicates the presence of a 
pool of relativistic electrons (CRe). Together with mag¬ 
netic fields this provides a favorable environment for high 
energy particle interactions between the CRes and the ICM 
which may be observable through the detection of 7 rays or 
hard X-rays [see, e.g. the excellent review by Brunetti and 
Jon es]|20l4l . However, due to the short diffusion times, 
CRes must be constantly replenished, e.g. through injec¬ 
tion by active galactic nuclei (AGN) or be created through 
secondary processes. 

Another intriguing possibility are hadronically-induced 
7 rays. Here CR protons (CRp) may be accelerated within 
the ICM through means of diffusive shock acceleration 
(DSA) and due to their large diffusion time remain within 
the cluster volume. CRp then interact with the ICM and 
produce 7 rays via decay of neutral pions. The latter has re¬ 
ceived particular attention as Pinzke and Pfrommer | 2010| 
have shown the emergence of a universal model describ¬ 
ing the CR interactions in a cosmological framework based 
on smooth-particle hydrodynamics simulations. In their 
model, the resultant 7 -ray spectrum is dominated by the 
aforementioned 7 r° decay and IC emission is essentially 
negligible. The spectrum shows the characteristic 7 r° bump 
at ^ 130 MeV and for energies > 500 MeV follows a 
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powerlaw with index 2.3 . 1 The resulting spatial distribu¬ 
tion is close to that of the thermal X-rays as it is expected 
that the CRs are following the gas. One key assumption 
when creating the model is that the maximum injection ef¬ 
ficiency, Cp,max at which protons can be accelerated via 
DSA is similar to that what has previously been observed 
in SNRs [Helder^f al. 2009]. 2 Together with the claimed 
universality of the spectrum, this allowed us to employ the 
joint likelihood technique, a statistical method in which 
each target is optimized according to its individual nui¬ 
sance parameters and then the information from each in¬ 
dividual likelihood is combined into a global likelihood by 
multiplying them [see Anderson 2014, this conference, for 
a technical discussion and applications]. 


The starting point for the study in Ackermann et al. 


| |2014a[ has been the extended HIFLUGCS catalog, a X 
ray flux-limited complete sample of nearby GCs. Selecting 
a set of 50 clusters, we found a global excess at the level 
of ^ 2 . 7 cr which however could be entirely attributed to 
previously non-detected individual cluster member galax¬ 
ies (with known counterparts in the radio band). Thus, 
with four years of LAT data, no cluster was detected and 
flux upper limits were set. The most constraining cluster 
in the sample is the Coma cluster with a reported integral 
flux limit of 4.0 x 10 - 10 ph/cm 2 /s assuming an extended 
emission characteristic according to the benchmark model 
by Pinzke and Pfrommer [20101 and 2.5 x 10 - 10 ph/cm 2 /s 
when considering a pointlike emission . 3 Based on the joint 
likelihood approach, we also find that in order to account 
for the non-observation, DSA must be either substantially 
less efficient (Cp,max < 21%) or conversely, the CR-to- 
thermal pressure ratio must be lower than 1 %, making the 
contribution of CRp’s to the ambient 7 -ray flux negligible 
[see also the discussion in |Vazza and Bruggen|2014| . 


3. Dark Matter constraints from Cluster 
Observations 


Given its non-detection, ongoing searches for 7 rays 
from GC can thus be used to constrain the available pa¬ 
rameter space of WIMP DM. Generically, the induced 7 - 
ray flux from WIMP pair annihilation can be expressed as 

MAn) = FJ^yL x f f p 2 ( r )didn'. (i) 

4?r JAn Jl.o.s. 

'-V-' S ---' 

$pp J—factor 


1 The interested reader is referred to ]Pinzke and Pfrommer||2010] l . The 
true spectrum is concavely shaped but for the considered LAT energies, 
it can approximated with a powerlaw as discussed in the main text. 

2 In their works, the authors adopt Cp,max = 50% as benchmark when 
calculating their 7 -ray predictions |Pinzke et a/.|2011|. 

3 These limits were calculated over the entire energy range of 

500 MeV to 200 GeV. 


In the above equation <f>pp refers to the particle physics 
term containing both the mass of the WIMP and its 
velocity-averaged pair annihilation cross section (cry). 
The second term is referred to as astrophysical, or J- 
factor and is the line of sight integral of the DM col¬ 
umn density. N-body simulations suggest that DM clus¬ 
ters (clumps) across all mass scales, forming sub haloes 
in addition to the smooth main halo. The amount of sub¬ 
structure as well as the properties are largely unknown and 
current N-body simulations do not have the capabilities yet 
to resolve the smallest substructures. Hence, extrapola¬ 
tions over several orders of magnitude are necessary. For 
this (abbreviated) discussion it is sufficient to address the 
amount of substructure by the introduction of a dimension¬ 
less boost factor b , which relates the J-factor obtained by 
assuming a universal NFW halo [Navarro et a/.||1997| to 
that obtained when considering different amounts of sub¬ 
structure. Depending on the extrapolation scheme, boost 
factors for clusters may vary between 0(30) to 0(1000) 
[see e.g. the discussion in Sanchez-Conde et al. ( [2011) . 
While the predicted annihilation flux profile is similar for 
both model predictions, the overall predicted flux may vary 
by orders of magnitude. This fact makes DM constraints 
from clusters (at least as far as annihilation is considered) 
more model dependent than e.g. those obtained from the 
observation of nea rby dwarf spheroidal galaxies |Acker- 
mann et a/.|2014b| . 4 


Ongoing work (based on a five year dataset) is focus¬ 
ing yet again on a subsample of the most massive nearby 
clusters, selected from the X-ray flux limited HIFLUCGS 
sample. We demand there to be no appreciable overlap by 
requiring a distance between each cluster of the sum of the 
virial radius of each cluster along with a 2 ° buffer account¬ 
ing for the tails in the LAT point spread function (PSF). 
For the resulting 34 clusters we construct spatial templates 
according to substructure models considering both a con¬ 
servative boost factor of 0(30) | |Sanchez-Conde and Prada| 
2014 ] and contrast this with a more optimistic choice of 
0(1000) [Gao et al 2012[ . For both configurations we 
perform a binned likelihood analysis. After having found 
the best fit parameters of our background fit, we construct a 
bin-by-bin likelihood function by assuming a simple pow¬ 
erlaw with index 2.0 in each energy bin which would ac¬ 
count for the cluster emission. The advantage is that the 
resulting flux limits can be used to directly test spectrally 
different models without the need of repeating the entire 
likelihood procedure [see Ackermann et al. 2014b for de¬ 
tails regarding the bin-by-bin method]. In Fig. [I] we show 
both the J-factor distribution from our chosen sample as 
well as the estimated sensitivity by selecting high-latitude 
regions in the sky which are selected such that the center 
a) does not contain a 3FGL source [The Fermi-LAT CoL| 


4 Note that the question of substructure is typically less important in 
the case of decay as the associated J-factor scales linearly with the DM 
density. 
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|laboration|2015[ and b) does not coincide with a cluster 
center or a circular region with a radius corresponding to 
the virial radius of the cluster . 5 


4. Challenging individual targets: very 
extended emission from the Virgo 
cluster 


While the discussed emission mechanisms vary appre¬ 
ciably with regards to the spectral form of the predicted 
emission, the studies that I discussed here have in common 
that the targets are large extended sources . 6 However, even 
among these extended sources, there are extreme cases. 
The largest target is the Virgo cluster, our closest neigh¬ 
bor which appears as a structure in the Northern part of 
the sky that extends up to 14° in diameter. The cluster 
itself is undergoing a complex merger in which the main 
clumps centered around the giant ellipticals M87 and M49 
are moving towards each other. 


5. Outlook 

Despite intense efforts no 7 -ray emission has been de¬ 
tected from clusters to date. In the meantime radio obser¬ 
vations have revealed more and more systems containing 
extended, yet faint radio sources. This radio emission will 
remain the driving force when searching for non-thermal 
emission from galaxy clusters. The LAT with its contin¬ 
uous sky survey capabilities will remain instrumental in 
probing the important ^ MeV — GeV domain which is 
much too low for air Cherenkov telescopes to be sensitive 
to. In particular the extension towards lower energies, en¬ 
abled by the latest reconstruction algorithms, collectively 
dubbed Pass 8 will help in achieving this goal [Atwood[ 
et al. es as it provides a better PSF and an increased 
acceptance towards lower energies. 

As for DM constraints, clusters will remain challenging 
targets - due to their large extension and the intrinsic uncer¬ 
tainties in their J-factors. However, they are also comple¬ 
mentary targets to probe in case evidence arises from more 
promising targets such as dwarf galaxies or the Galactic 
center. 
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to emphasize that the model for the Galactic foreground 
emission that is usually employed when analyzing Fermi- 
LAT data is optimized for point source searches. Indeed, 
when confronting a dataset comprising three years of col¬ 
lected photons between 100 MeV and 100 GeV, we find an 
extended excess if we em ploy the standard diffuse model 
jAckermann et a/.|2015| . 7 However, when systematically 
performing a scan over the entire ROI by using a uniform 
disk of 3° radius, we find two distinct maxima which are 
spread out and appear away from both sub clump cen¬ 
ters as shown in Fig. [2] Moreover, when using a set of 
alternatively derived diffuse foreground models [ [Acker- 1 
|mann et al. |2012 de Palma et al. |2013) , the significance of 
this extended excess varies appreciably, implying that the 
source of the excess may be due to an incomplete modeling 
of the Galactic foreground emission. 


5 For this analysis we select the subset of P7REP photons that pass the 
CLEAN class and apply the recommended instrument response function 

P7REP_CLEAN_V15. 

6 With large we refer to an emission radius of ~ 2 — 3° (as in the case 
of the Fornax and Coma cluster, respectively). 

7 For this analysis we select Pass 7 (P7V6) photons passing the 
SOURCE selection together and apply the recommended models for 
modeling the Galactic diffuse and isotropic emission. The reader is 
referred to the web pages of the Fermi Science Support Center for 
details: http: / / f ermi .gsfc.nasa. gov/s sc/data/access/ 
lat/BackgroundModels.html 
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Figure 1: Left: Expected sensitivity (evaluated by calculating the combined one-sided 95% C.L. upper limit) for annihilating WIMP 
DM into bb from a combined analysis of 34 galaxy clusters in 5 years of LAT data using photons from 500 MeV-500 GeV. Each 
cluster model contains a conservative amount of substructure. The red solid line corresponds to the first LAT team publication based 
on 11 months of data [ Ackermann et al. 2010] in which each individual limit has been weighted with the assumed J-factor (NFW 
only). Right: Distribution of J-factors for our chosen cluster sample. In the main text we refer to the setup by|Sanchez-Conde and 
|Pmda|pOT4| as conservative setup. The optimistic configuration refers to the distribution labeled with |Gao et al. |2012|. 



Figure 2: Grid scan of TS values when placing a uniform disk 
emission template with a powerlaw and 7 = 2.0 on various 
positions in the map. This representation allows to map the 
observed excess when the foreground Galactic diffuse emission 
is modeled using the standard template for point source analysis. 
For reference we add the reported NED positions for the central 
galaxies M87 and M49 as diamond and star-shaped markers, 
respectively. 
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